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Abstract. These lectures are devoted to elementary particle physicists and assume 
the reader has very little or no knowledge of cosmology and astrophysics. After a brief 
' historical introduction to the development of modern cosmology and astro-particles 

in which the Hot Big Bang model is defined, the Robertson- Walker metric and the 
dynamics of the Friedmann-Robertson- Walker cosmology are discussed in section 2. 
In section 3 the main observational features of the Universe are reviewed, including a 
| description of our neighbourhood, homogeneity and isotropy, the cosmic background 

radiation, the expansion, the age and the matter content of the Universe. A brief 
account of the thermal history of the Universe follows in section 4, and relic abundances 
are discussed in section 5. Section 6 is devoted to primordial nucleosynthesis, section 
7 to structure formation in the Universe and section 8 to the possibility of detection of 
. the dark matter in the halo of our galaxy. In the relevant sections recent developments 

are included, such as several so called "crisis" (the age crisis, the cluster baryon crisis 
and the nucleosynthesis crisis), and the MACHO events that may constitute the first 
detection of dark matter in the halo of our galaxy. 
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^ . 1. Historical Introduction 

Cosmology, from the greek noafiOi;, that means "order" , is the study of the origin 
and evolution of the Universe. The initial point of modern cosmology may be taken 
to be 1905 with the advent of Einstein's Relativity Theory, extended to General 
Relativity in 1915. The first solutions to Einstein's equation for a homogeneous and 
isotropic but time evolving Universe were found by Friedman in 1929 and in 1935- 
1936 Robertson and Walker showed that Friedman's metric could be derived from 
homogeneity and isotropy alone, independently of General Relativity. Friedman's 
equations show that a Universe containing matter and radiation is expanding (or 
contracting). In 1929 Hubble found the first indication of the recession of faraway 
galaxies with a velocity v proportional to its distance d, v = H a d, where H is 
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called the Hubble constant, the present value of the Hubble parameter (constant 
in space but not in time). Notice that the subscript zero indicates a present value 
in cosmology. 

Together with the belief that we are not in a special place in the Universe, elevated 
to a principle in modern cosmology, the Cosmological Principle, the recession from 
us of distant galaxies means that the Universe is expanding (so that any two bodies 
not gravitationally bound to each other are receding from each other). If the 
Universe was smaller in the past and it expanded adiabaticlly, it was denser and 
hotter. Gamow was first in taking seriously this idea. He proposed in 1946 that 
all the matter in the Universe was hot and reached nuclear and larger densities in 
the past. This is the "Hot Big Bang" model [1]. This model implies a beginning of 
the Universe a time t Q ~ H^ 1 ago. Hubble's estimate of H , the Hubble constant, 
was around 5000 km/(sec Mpc), larger than the presently known value of 40-100 
km/(sec Mpc), what yielded a lifetime of the Universe t ~ H~ l ~ 4.5 x 10 9 yr, 
shorter than the already known age of the Earth (1 pc = 1 parsec = 3.2615 light-yr). 

In part because the measured t Q was too short at the time, in part due to the 
philosophical difficulty of accepting that the Universe had a beginning, the com- 
peting "Steady-State Universe" model was introduced by Bond, Gold and Hoyle 
in 1948. In this model the Universe is homogeneous also in time, it always looked 
as it looks now. Because in the Hot Big Bang model matter and radiation were 
in equilibrium at early times, a remanent of radiation with a black body spectrum 
is expected (that would not be there in a Steady-State Universe), as predicted by 
Alpher and Herman [2] in 1949. The observation of this radiation by Penzias and 
Wilson [3] (as explained by Dicke et al. [4]) in 1965, confirmed the Hot Big Bang 
model. It is remarkable that Alpher and Herman predicted the temperature of this 
cosmic microwave background radiation, CMBR, to be 5°K, quite close to its ac- 
tual temperature of 2.7°K. Actually, the experimental confirmation of the Hot Big 
Bang already existed when Gamow first proposed it, in 1946. Astronomers knew 
since 1941 that interstellar CN molecules are excited in high rotational levels, as 
seen through CN absorption lines in the spectrum of stars, indicating that those 
molecules are in a thermal bath at about 3°K. 

Gamow's main interest in proposing the Hot Big Bang was the primordial syn- 
thesis of nuclei. It was soon realized that the synthesis of heavy elements had to 
be done in stars, what lead to the study of stellar evolution. Only in 1966 the 
quantitative study of primordial nucleosynthesis started, the synthesis of 4 He at 
first and that of D, 3 He and 7 Li only later, in the early 70's. 

This period, the early 70's, can be considered as the beginning of astro-particles. 
In 1967 (the year when the standard electroweak model was proposed) Sakharov 
noted the essential elements necessary to generate the asymmetry between mat- 
ter and antimatter in our Universe, namely the generation of baryon number or 
baryogenesis. One of them is the violation of baryon number, that encountered in 
Grand Unified Theories (first proposed by Georgi and Glashow in 1974) its most 
natural source. In the 70's cosmological and astrophysical implications of neutrinos 
were first studied, leading to the battery of tests that limit the number of neutrino 



species, their mass, cosmological density, lifetime, decay modes etc. Also, it was in 
1975 when the first evidences for the dark matter were discussed and the idea of 
inflation came about in 1980. (For references and more details about the content 
of this section so far, see e.g. [5]) 

The last fifteen years can be considered to be the golden age of astro-particles, in 
which the connection of particle physics and cosmology has lead to a number of im- 
portant ideas. Among them is the idea of non-baryonic dark matter -see section 8-, 
that is explored in simulations of galaxy formation confronted with observational 
cosmological data and in experimental dark matters searches, both with cryogenic 
detectors and large neutrino detectors, supplemented by accelerator bounds. Other 
examples are the models for baryogenesis and for inflation, the interplay between 
the properties of proposed particles and their cosmological and astrophysical con- 
sequences, the introduction of particle physics experimental techniques to study 
astrophysical problems (such as in large "neutrino telescopes" or in the search for 
MACHOS -see section 8- in the halo of our galaxy). As an example of this last 
point, let us mention that experiments designed to search for the decay of protons 
predicted by Grand Unified Theories (1MB and Kamiokande) observed for the first 
time neutrinos from a supernova, SN 1987A, inaugurating with this observation of 
neutrinos from outside the solar system, what we can call neutrino-astronomy. 



2.1 The Model - The Hot Big Bang (Hot BB), the standard model of cosmol- 
ogy, establishes that the Universe is homogeneous, isotropic (thus, it is a Friedmann- 
Robertson- Walker Universe) and expanding from a state of extremely high tem- 
perature and density. Thus the early Universe can be described as an adiabaticlly 
expanding classical gas of relativistic particles, namely radiation, in local thermal 
equilibrium at a temperature T, that changes with time, T(t). The lifetime of the 
Universe t is counted from the moment the expansion started, taken to be t — 0. 

This model is based on General Relativity, the Cosmological Principle and three 
major empirical pieces of evidence, namely the Hubble expansion, the cosmic black- 
body microwave background radiation (CMBR) and the relative cosmic abundance 
of the light elements (up to 7 Li). Because the model is based on General Relativity, 
it is for sure not valid in the realm of Quantum Gravity, T > M P i ank = 1.22 x 10 19 
GeV and t < 10~ 43 sec. The Cosmological Principle postulates that we do not live 
in a special place in the Universe, by requiring that every comoving observer in the 
"cosmic fluid" has the same history. The "cosmic fluid" has as particles clusters of 
galaxies, and "comoving" in practice means at rest with the galaxies within a 100 
Mpc radius (one parsec, 1 pc = 3.26 light-years = 3.09 x 10 18 cm), as will become 
clear below (see section 3.2). The Hubble parameter H provides the proportionality 
between the velocity v of recession of faraway objects, and their relative distance 



2. The Hot Big Bang Model 
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The CMBR was produced at t rec ~ 3 x 10 5 y, the recombination epoch, when atoms 
became stable. It has a blackbody spectrum and it is remarkably isotropic. Finally, 
the earliest available proof of the consistency of the Hot BB model is provided by 
the cosmological abundance of 4 He and of the trace elements D, 3 He and 7 Li. Their 
abundances, differing by several orders of magnitude, are well accounted for in 
terms of nuclear reactions that occur at £ B bn — 10~ 2 — 10 2 sec, the Big Bang 
nucleosynthesis (BBN) epoch, when the photon temperature was Tbbn — 10 — 0.1 
MeV (necessarily below the binding energy of the light nuclei). Let us expand 
on these points. 

2.2 Friedman- Roberts on- Walker Models Homogeneity and isotropy re- 
strict the space-time metric to be of the Friedmann- Robert son- Walker form, 



This metric depends on a global scale factor a(t) that describes the overall expan- 
sion of the Universe, i.e. a physical distance A increases with the Hubble expansion 
as A = a(t)X c , where A c is the constant distance measured in a comoving refer- 
ence frame (a frame that follows the Hubble expansion, like a grid painted on an 
expanding balloon), whose coordinates are r, 9 and 0. So \ c ,r,9,(f> do not change 
with the expansion of the Universe. With the usual choice of a Q —l for the present 
value of the scale factor, A c corresponds to the present physical distance. The time 
coordinate t in Eq. (2) is the proper time measured by a comoving observer, an 
observer (located at fixed r, 9 and (f>) that moves with the local Hubble flow, k is 
the curvature parameter. The spatial curvature, 3 K(t) = k/a(t) 2 , can be either 
positive (corresponding to a closed Universe with the geometry of a 3-sphere), or 
zero (corresponding to a flat open Universe), or negative (corresponding to an open 
Universe with the geometry of a 3-hyperboloid) depending on the value of the cur- 
vature factor. Through rescaling of the coordinates, k can be chosen to be +1, 
or —1 respectively. 

This metric yields the Hubble expansion if the Hubble parameter H(t), defined 
as H(t) = a(t)/a(t), is positive, since taking t, 9, and constant, ds = dD(t) where 
D(t) is the proper distance, D(t) = a(t) f dr(l — A;r 2 )~ 1//2 , so that the proper 
velocity V = dD(t)/dt follows Hubble's law, V = HD. However V and D are 
not operational quantities (those that can be measured) such as the velocity of 
recession v and the luminosity distance di (see, for example, [5]). v is measured 
through the red-shift z of the observed light (of wavelength A D ) with respect of the 
emitted light (of wavelength A e ), 1 + z = A D /A e = a(t )/a(t e ) ~ v/c, where c is 
the velocity of light. This is the non-relativistic expression, only valid for 2 < 1, 
otherwise (v/c) = [(1 + z) 2 — 1] [(1 + z) 2 + l] -1 . g?l is measured through the absolute 
luminosity C of the emitting body (whose difficult determination is the origin of 
the uncertainties in the measurements of H) and the measured flux J r ) d\ = C/ 
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reduces to Eq.(l) with d = d L) for v <C c. The deceleration parameter q Q = 
— (aa/d 2 ) , if not zero, becomes important at moderate z, where yet poorly known 
changes in the sources due to their evolution also become important. For this 
reason q Q is only known to be at most of order one. H a , the Hubble constant, the 
present value of the Hubble parameter H, is at present observationally determined 
within a factor of 2 to be 

km . , 

H = 100 h z-— , 4 

sec Mpc ' v ; 

with h— 0.4 -1. Through Hubble's equation, recession velocities or red-shifts are 
translated into distances, 

d L ~ 10~ 2 ; — h^Mpc = 3000 z h^Mpc , (5) 

km/sec 

for z < 1. Notice that the units of distance are (/i _1 Mpc) in this case. 

2.3 Evolution of the Universe - The dynamics of the global scale factor 
a(t) is determined by the content of the Universe (in matter, radiation and vacuum 
energy) through the Einstein field equations, R^ v = \Rg ilv = 8nGT^ u — Ag^. 
Spatial homogeneity and isotropy require the total stress-energy tensor to be 
diagonal and with equal space components. The simplest realization corresponds 
to a perfect fluid = diag (p(t) , —p(t) , —p(t) , —p(t)) where p(t) and p(t) are 
the energy density and the pressure respectively. A, the cosmological constant 
introduced by Einstein in 1917 in order to obtain a static Universe (idea rejected 
experimentally by the discovery of the Hubble expansion in 1929), is associated 
in modern terms with vacuum energy. It can be incorporated into by adding 
Pvac = A/87TG and p vac = —p VSLC to the p and p of radiation and matter to obtain 
Ptotai = P + Pvac and ptotai = P + Pvac- With this there are only two independent 
Einstein equations, from the 00 and ii components respectively. From p — 0, v — 
we obtain the Friedmann equation, 
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whose Newtonian meaning of energy conservation for a unit mass test particle can 
be seen (for A = 0) by writing it as d 2 /2 - GM/a = -k/2 for M = (4/3) na 3 p, 
where the total energy is —k/2. Combining the 00 and ii components one gets 

d -AttG , . A 

- = ^^(P + 3p) + -, (7) 
a 3 3 

that for matter (gas, dust), for which p = (and A = 0), has the meaning of 
acceleration equal force over mass, a = —GM/a 2 . 

In order to determine the three unknown function of time in these equations, a(t), 
Ptotai (t) and Ptotai(0 we need a third independent equation. This is the equation 



of state, p = p(p) provided by thermodynamics. For only relativistic particles 
and radiation, pr = prc 2 /3, while for only matter (gas, dust) pu = (for a gas 
Pm — n kT ~ Pmv 2 *C pc 2 , where n is the number density and k is the Boltzmann 
constant). It is convenient to introduce the energy conservation equation, T^ v = 
for /i — 0, that is not independent of Eqs. (6) and (7), p = — 3(d/a)(p +p). Using 
the equations of state in this equation, one obtains p as a function of a, pr ~ a~ 4 
and pm ~ cl~ 3 , and using these in Eq. (6) one obtains, 

a R ~ t 1 / 2 , = 1 a M ~ t 2/3 , E M = \ t - (8) 

Thus, given a certain value of H a , a Universe radiation dominated for most of its 
life is younger than a matter dominated one. The relations in Eq. (8) can be easily 
obtained by using a ~ T~ 4 (this is the entropy conservation condition, derived 
from the energy conservation equation in thermal equilibrium). Then, we know 
that t^ 1 ~ H ~ yfp (from the Friedmann Eq. (6)) and we know that Pr ~ T 4 and 
p M ~ m/a 3 ~ T 3 , therefore t ~ T~ 2 ~ a 2 for radiation and t ~ T~ 3/2 ~ a 3/2 . 

With only vacuum, namely with p = p = and A > constant, neglecting the 
curvature term we obtain from Eq. (6) an exponential expansion of the scale factor 
a ~ e HT with H = A/^/3. Notice that with an exponential expansion the curvature 
term becomes negligibly very fast, k/a 2 A/3. This is precisely why a period of 
vacuum dominated expansion in the early Universe, called "inflation", has been 
proposed to explain the flatness (see section 3.7) of the Universe. 

Inflation would happen during a period of supercooling of the Universe (so that 
the density of radiation and matter, proportional to T a , a = 4, 3 become negligible) 
while the Universe is rolling down the almost constant potential V{4>) ~ V Q of a 
scalar field, the "inflaton" 0, so that p vac ~ V Q ^> Pr,Pm- At the end of inflation 
it is necessary to convert V Q into radiation and matter, through the decay of the 
4> field, in a process called "reheating", that generates the Universe as we know it, 
still with a negligible curvature. We can rewrite Eq. (6) in the form 

k Ptotal -y Ptotal -y ^ j ^g-j 



H 2 a 2 (3H 2 /8ttG) " p c 

where Q is the density in units of the critical density, p c = 3H 2 / (871G) = 
10.5 h 2 (keV/cm 3 ) = 1.88 l(T 29 /i 2 (#/cm 3 ). Notice that spatially open, flat and 
closed Universes (with k < 0, k = and k > respectively) correspond to Q < 1, 
Q = 1 and Q > 1 respectively. So, a long enough period of inflation produces 
a~ 2 — > 0, thus 1 — Q^ 1 ~ (k / a 2 ptotal) — (since ptotal = Pvac is constant). That is, 
n — 1 after a long enough inflation. 

Notice that only for A = there is a correspondence between spatial geometry 
and the future evolution of the Universe, at t — > 00. Let us call here fi G the density 
of matter and radiation, Q Q = Q R + Q M , so that the total density is Q = Q Q + JIa- 
Because for matter and radiation p ~ a~ a with a = 3 or 4, for large values of the 
scale factor a the curvature factor k/a 2 dominates the r.h.s of Eq. (6). In this 



case, an open or flat Universe, i.e. with Q = fi < 1 (k < 0), expands forever and 
a closed Universe, with il = Vl > 1 (A; > 0) recollapses. A cosmological constant 
A dominates instead the r.h.s. of Eq. (6) for large values of a, and for A > the 
Universe expands exponentially forever at large times (even if Q > 1), as can be 
seen in Fig. 1 (taken from [6]). 

Notice that for matter or radiation p + 3p > 0, what, through the Eq. (7) (with 
A = 0), implies that a < 0, so the expansion of the Universe is slowed down by 
the gravitational pulling of the mass content of the Universe, and consequently 
t < H^ 1 (since t = H^ 1 corresponds to a — 0). On the other hand, for a vacuum 
dominated Universe the pressure is negative, p v3uC = — p vac and Eq. (7) shows that 
gravity is repulsive, a = Aa/3 > and consequently t > if -1 . Thus, a non-zero 
cosmological constant yields larger values of t Q for a given H Q (see Fig. 1) than 
matter and radiation alone. We see, therefore, how a non-zero A can resolve the 
present possible conflict between a large value of H Q and a large value of t Q , the 
"age crisis" we will mention later (see section 3.6). (For references and more details 
about this section see e.g. [5]) 



3. Observational Features Of The Universe 

3. 1 Topography of Our Neighbourhood - Galaxies are the building blocks 
of our Universe. Normal galaxies have masses of 10 8 to 10 12 M & (where M & is the 
solar mass), while dwarf galaxies have 10 6 to 1O 7 M (same mass as that of globular 
clusters, old systems of stars spherically distributed around the center of galaxies 
like ours). We live in a spiral galaxy, the Milky Way, with a thin disk of a few 100 



FIGURE 1. Evolution of the cosmic scale factor a(t) in different FRW models. 



pc of thickness and approximately 12 kpc of radius. The Sun is in one of its arms, 
at about 8 kpc from the center. 

Galaxies form gravitationally bound systems: binary systems, groups (that are 
systems of a few to 10 galaxies, with a typical size of 1 Mpc), clusters (systems 
with hundreds of galaxies extending from a few to 10 Mpc) and superclusters (with 
thousands of galaxies in a radius of 10-50 Mpc). We live in a group, the Local 
Group, together with another large spiral galaxy, Andromeda, at roughly 700 kpc 
of the Milky Way, and several small galaxies. 

At the scale of superclusters, i.e. scales of 20-100 Mpc, the Universe is better de- 
scribed in terms of "walls", "filaments" and "voids" [7]. There are few galaxy maps 
at even larger scales, i.e. scales larger than 100 Mpc, and they show smoothness 
at those scales. Present maps only cover a small fraction of our visible Universe, 
whose radius is ct D ~ 10 10 light yr ~ 3000 Mpc. 

While Hubble's law and red-shift measurements are the only tools to obtain the 
distance of very faraway objects, the distance to nearer ones can be estimated 
separately from their red-shift, what allows to measure "peculiar" velocities v pec , 
i.e. velocity components due to gravitational acceleration that add to the local 
velocity due to the Hubble flow H Q d to give the total velocity v = v pcc + H d. 
In this way it has been determined that large regions of about 100 Mpc of our 
local neighbourhood are moving towards a very massive system, called the Great 
Attractor, with large peculiar velocities of about 600 km/sec [8] [9]. 

3.2 Homogeneity and Isotropy - Although these are assumptions incorpo- 
rated into the BB model by the Cosmological Principle, there are strong observa- 
tional foundations of the homogeneity and isotropy of the Universe at scales larger 
than 100 Mpc. 

As we have just seen, the Universe looks lumpy and large regions move with 
large peculiar velocities at scales A < 100 Mpc, however the Universe looks smooth 
at larger scales and we also know it was smooth in the past. There are few 3- 
dimensional galaxy maps that reach the necessary large scales to test homogeneity, 
containing in total only of the order of 10 4 galaxies, but by the end of the century 
10 6 will be mapped. These are red-shift surveys where distances are obtained from 
z though Hubble's law. Angular photometric galaxy surveys instead, provide 2- 
dimensional information for several 10 6 galaxies. Both types of galaxy maps show 
smoothness at A > 100 Mpc, in the case of 2-dimensional maps by showing that 
angular correlations decrease steeply with increasing angular distance in the sky 
(see for example [6] for references). But the best evidence of homogeneity and 
isotropy is provided by the isotropy of the CMBR, for which measurements of the 
temperature in different regions of the sky show very small anisotropies, smaller 
than 10 -4 , i.e. the r.m.s. temperature fluctuation measured between two points 
in the sky separated by an angle 6 is (ST/T)g ^ 10~ 5 [10] [11]. This shows both 
the smoothness of the Universe at the moment of emission of the radiation and the 
isotropy of the expansion since then, because the present measured wavelength of 
the CMBR photons now is A now = a(t)X e , where A c is the wavelength at emission, 
and any space dependence of the scale factor a(t) would produce anisotropies. The 



CMBR photons are emitted from the last scattering surface, that is the Universe 
at recombination (i~3x 10 5 yr, T ~ 0.3 eV and z ~ 1100). 

3.3. Cosmic Microwave Background Radiation - The best blackbody in 
the Universe is the Universe. The COBE (Cosmic Background Explorer) Satellite 
measured the CMBR spectrum in 1992 for wavelengths between 0.05 and 1 cm, 
and found a blackbody with temperature T Q = 2.726 ± 0.010°K, with deviations 
of less than 0.03% [12]. This provides fundamental evidence for a Hot BB. As we 
just mentioned, the CMBR was produced at the recombination epoch, when atoms 
became stable and replaced ions and electrons in a plasma as the constituents of 
the Universe. At this moment the mean free path of photons, that was very short 
in the preceding plasma, became very long, so that the photons that last scattered 
then can reach us. Knowing so well the CMBR temperature, we know with great 
accuracy the number and energy density of the CMBR photons, that are the most 
abundant in the Universe by several orders of magnitude (see e.g. [5] p. 143), 
n 7 = 2C(3)T 3 /vr 2 = 411/cm 3 , p 7 = tt 2 T 4 /15 = 4.71 x 10- u (g/cm 3 ) Using p c in Eq. 
(9), we see that fi 7 ~ 10 -5 . 

Anisotropies in the CMBR are expected due to our motion with respect to the 
CMBR rest frame, and due to the density inhomogeneities that triggered structure 
formation in the Universe. In fact COBE measured a dipole anisotropy correspond- 
ing to a velocity of 627 ± 22 km/sec of our Local Group of galaxies with respect 
to the CMBR rest frame (COBE even saw the rotation of the Earth around the 
Sun!), and measured anisotropies (5T/T)g for angles 9 = 7° to 90°. At 90° COBE 
measured a quadrupole anisotropy of (5T/T) 90 ° — 0.5 x 10 -5 [10]. Results from 
other experiments in balloons are available now at smaller angles, 9 = 0.5° to 90° 
and the results show 8T/T < 10 -5 after subtracting the dipole. 

The horizon size at recombination, ct rec , corresponds to an angle 
9 H ~ l°(r2) -1 / 2 in the present sky. Larger angles correspond to causally dis- 
connected regions at the time of emission of the CMBR photons. Thus the 
smallness of the anisotropies at angles 9 < 1° indicates that the Universe was 
smooth at recombination, and at larger angles, 9 ^ 1° it indicates that the Uni- 
verse is smooth on very large scales. The temperature autocorrelation function, 



C(9) = < (5T(a)/T) (5T(a + 9)/T) > a , where 5T/T = (T-T)/T are 



temperature fluctuations with respect to the average temperature T = 2.726°K, 
is computed by measuring 5T/T at some position a in the sky (a is given in 
spherical harmonics by two angles), multiplying it with 5T/T at another position 
(a + 9) separated from the first one by an angle 9 and averaging over all positions. 
The expansion of C{9) in Legendre polynomials defines the multiple moments Cg, 
C(9) = (47r) _1 X^(2£ + l)QP^(cos#), so that the anisotropy at a certain angle 9 is 



We see that £ > 200 correspond to scales within the horizon at decoupling, where 
matter can move due to density perturbations and where the photons pick up 





(10) 



additional energy due to the scattering from moving matter (mainly Thompson 
scattering on electrons). This gives origin to Doppler (or acoustic) peaks in the 
spectrum of multiple moments. The position of the first peak should happen at 
the scale of the horizon £ pea k — 200(f2) _1 , as shown in Fig. 2 [13] for different dark 
matter models (see section 7). Thus, in the future, the position of the first peak 
will determine f2 within a 10%. The height of the peaks depends on the density 
in baryons, Qb and the location and relative height of the peaks will allow to dis- 
criminate different models of galaxy formation, once all Cg are measured. Actually 
only the Cg for t < O(10 3 ) are relevant, because larger values of £ correspond to 
scales 9 < 8° that are inside the thickness of the last scattering surface, so the 
information at those small scales is smeared out [11]. 

3-4 The Expansion of the Universe - Hubble's law in Eq.(l) is valid at 
z <C 1. At z 1 the law becomes quadratic (see Eq. (3)) and the deceleration 
parameter q Q measures the rate at which the gravitational attraction is slowing 
down the expansion. Actually q a = Q Q /2 — ft\ (for matter dominated Q ), so for a 
large cosmo logical constant component fl\, q Q could be negative and the expansion 
could be accelerated (recall that gravitation is repulsive in a vacuum dominated 
Universe). Not much can be said observationally about q a yet. So let us concentrate 
on the Hubble constant H Q . 

Redshifts can be measured accurately. This allows to determine the recession 
velocities v. In general they have two components v = H Q d + v pec . Thus recession 
velocities can be attributed solely to the Hubble flow at distances large enough 
for the peculiar radial velocities v pec , the velocity components due to gravitational 



FIGURE 2. CMBR anisotropy power spectra predicted by four models (lines) of structure 
formation (see section 7). The band shows the expected experimental error in future satellite 
experiments. 



acceleration, to be negligible, H Q d ~ v ^> v pec , what happens at d ^> /i" 1 Mpc, say 
distances larger than 50 Mpc, because t> pcc ~ 0(100 km/sec). 

In order to measure H Q = v/d (Eq. (1)), distances and velocities must be deter- 
mined separately. The dificulty in the measurement of H Q reside in determining d. 
Different techniques to measure distances lead to the spread of roughly a factor of 
two in results for h (the Hubble constant in units of 100 (km/sec Mpc) see Eq. (4)). 
Most methods give large values of h, h ~ 0.8. For example, a recently publicized 
result of the Hubble Space Telescope is h — 0.8 ± 0.17 [14]. This determination is 
based on the observation of twenty Cepheids in the Virgo cluster. Cepheids are vari- 
able stars, whose period of variability is associated with their intrinsic luminosity, 
thus, the apparent luminosity tells the distance to the star. Some determinations 
give smaller values of h. For example, using type la supernovae a determination 
of h = 0.5 ±0.1 was made (see e.g. [9] for references). 

The problem with a large value of h, h >^ 0.65, is that it yields a too young 
Universe, unless A ^ or Q Q < 1, namely unless we live in a Universe with positive 
cosmo logical constant, or open, or a combination of the two. This goes against 
aesthetic beliefs in a flat, matter dominated Universe, giving origin to what some 
call an "age crisis" (see section 3.6 below). 

3.5 The Age of the Universe - There are three main methods to deter- 
mine t . The technique of nuclear cosmochronology, uses the relative abundance of 
radioactive isotopes at present and at production (through rapid neutron capture 
processes, r-processes, most probably in supernovae). Comparing the two gives 
t = (10 - 20) x 10 9 yr [15]. 

White dwarfs are the last stage of low-mass stars. Since they are faint they can 
only be seen near the solar system. There is a drop in the number of white dwarfs 
with luminosity smaller than 3 x 10~ 5 of the solar luminosity. Combining cooling 
models with the assumption that this cut indicates that there are no older white 
dwarfs that had the time to become fainter, gives the age of the disk of our galaxy. 
This method gives t a = (9 - 10) x 10 9 yr [16]. 

Globular clusters are old systems of a few-million stars that formed all at the 
same time. The evolution of stars depends on their initial mass, the more massive 
have shorter lifes, and finish burning Hydrogen earlier. The age of globular cluster 
is determined by measuring the absolute luminosity of the stars that are finishing 
burning H, using well known stellar evolution theory. The lifetimes obtained in this 
way are t Q ~ (13 — 15) x 10 9 yr and it is very difficult to get to t Q < 11 — 12 x 10 9 
yr [17]. Thus, t Q >^ 13 Gyr is taken at present to be a reasonable lower bound of 
the age of the Universe and t a >^ 10 Gyr is taken to be an absolute lower bound, 
a minimum age already uncomfortable to accommodate globular clusters. 

3.6 Age Crisis? — The parameters t , H Q and f2 are not independent, actually 
H t = 1.06 (/i/0.80)(t o /13 x 10 9 yr) = ft(Si), where ft is a function of Q that 
depends on the content of the Universe (see e.g. [5]). For a flat f2 G = 1 Universe 
matter dominated by the dark matter (see section 3.7) with A = (until recently the 
model preferred by most cosmologists due to its simplicity and aesthetic appeal), 
one has H t = 2/3 what implies that t >^ 13 Gyr requires h < 0.50 and t >^ 10 



Gyr requires h < 0.65. As mentioned above, in section 3.4, these values of H Q are 
lower than most present determinations (but not all). If h is actually larger than 
0.65, then we live in a Universe with a non-zero cosmological constant, or open, or 
both. 

If our Universe is spatially flat, namely Q = Q Q + Qa = 1, the frequency of 
gravitational lensing of quasars by nearer galaxies gives a bound of £l\ < 0.7 [18]. 
The number of lensings depends on A, because with A > the distance to a quasar 
of a given redshift is larger than for A = 0. With Q = 1 and < 0.7 one gets 
H t < 0.96. In this case t > 13 (10) Gyr requires h < 0.72 (0.94). On the other 
hand if A = and h > 0.65, we need Q(— Q Q ) < 1, namely we live in an open 
Universe. For a matter dominated Universe H t increases for decreasing fl Q , from 
2/3 for Q Q = 1 to 1 for Q D = (namely for an empty Universe t = H' 1 , since 
a — 0). Dark matter measurements indicate that there is at least Q Q = 0.1 in 
matter, what means H t < 0.9. In this case (with A = 0) t Q > 13(10) Gyr requires 
h < 0.68 (0.88). Of course a combination of A ^ with VL < 1 would also be 
possible, but even less appealing. 

Why does A ^ or Q < 1 seem so unappealing? There is a problem of fine 
tuning in both cases. A cosmological constant corresponds to a non-zero vacuum 
energy. We cannot explain why the vacuum energy should be zero at present, but 
it is even less appealing to explain why if non-zero it is so small with respect to 
characteristic values of the potential energy of scalar fields in particle physics. Take 
as an example the usual Higgs potential responsible for the spontaneous breaking 
of the electroweak symmetry. Take to be the usual Higgs doublet, then the 
potential energy density is V = \(\(j)\ 2 — v 2 ) 2 and the hight of this potential is 
V = V(<j> = 0) = Xv 4 = A(100 GeV) 4 . In units of GeV 4 the critical density is 
p c = lO" 46 GeV 4 , so Pvac = p A < Pc = lO-^A-Vj 

With Q ^ 1 there is an equally important fine tuning to take into account. If 
Q 7^ 1, Q increasingly separates from 1 as t increases in a Universe radiation or 
matter dominated. Dividing Eq. (9) by Q, we find 1 — = 3i^(87rGa 2 p) _1 = x(t) 
and x(t) ~ a for matter ~ a 2 for radiation (see Eq.(8)). Thus, in order to have 
Q = 0(1) (but ^ 1), n should have differed from 1 in the past by an extremely 
small amount, at most x(t) — 10 -10 at the nucleosynthesis epoch (T ~ 1 MeV) 
or 10~ 32 at the Planck time (T ~ 10 19 GeV). These small initial values of x(t) are 
unappealing. 

3. 7 Matter Content of the Universe - The luminous mass, namely the 
matter associated with typical stellar populations, accounts for Lum ~ 0.003/;, -1 
[19]. It is determined through luminosity measurements knowing the average mass 
over luminosity ratio. This is, for example, 5 times the solar ratio in the solar 
neighborhood and 12h times the solar ratio in the core of elliptical galaxies. 

However there is much more than the luminous mass in the Universe. It is by 
now a well established fact that on all scales larger than about 10 kpc there is 
a discrepancy between the amount and distribution of the luminous mass, and 
the mass detected through its gravitational effects. The gravitationally dominant 
mass component of the Universe is "dark", i.e. it is not seen either in emission 



or absorption of any type of electromagnetic radiation. This is called dark matter 
(DM). The most robust evidence for the DM is given by the rotation curves v(r) 
of spiral galaxies, i.e. the orbital velocities of stars and gas clouds orbiting a spiral 
galaxy at a distance r from its center. These velocities remain constant outside 
the region where the light falls exponentially off. Since v(r) 2 = GM(r)/r, the 
flatness of the rotation curves v(r) beyond a few kpc from the center means that 
the mass contained within a radius r grows linearly with r, M(r) ~ r. This 
indicates the existence of a quasi-spherical DM halo around each galaxy whose 
density falls off as r -2 outside a core radius of a few kpc. These DM haloes amount 
to Vl DM = 0.02 — 0.05 [20] [19]. The first indications of the necessity of DM came 
from the study of the Coma cluster by Zwicky in 1930. Applying the virial theorem 
to clusters of galaxies the total amount of matter in the cluster can be inferred from 
the characteristic velocity and average separation of its galaxies. The DM measured 
in this manner is Odm — 0.1 - 0.2 [19]. A more recent method uses the gravitational 
lensing of very far galaxies by clusters, whose total mass can be reconstructed from 
the deformation of the lensed galaxies into arcs and arclets. This method also yields 
f2 DM ~ 0.1 - 0.2 [21]. Peculiar velocities are due to gravitational accelerations due 
to masses nearby. Large scale peculiar velocity flows, at scales of 50 Mpc and 
larger, can then be used to infer total mass distributions, what yields estimates in 
the range Odm ^0.2-1 [9]. 

With respect to the density in baryons (namely nuclei) Qb, Big Bang nucleosyn- 
thesis (BBNS) predicts the density of light elements in agreement with observations 
only for Vtsh 2 ~ 0.01 — 0.02 (see more about nucleosynthesis in section 6), what 
for h = 0.4 — 1 means 0.01 < Qb ~ 0.13. This means that all the DM in the halo 
of galaxies could be baryonic (in the form of macroscopic objects, such as failed 
stars or primordial black holes for example), but also may indicate that some of 
the baryons are dark (Qb > ^Lum) if h < 1. 

The BBNS estimate of tt B combined with the measurement of a "large" amount 
of gas in rich clusters of galaxies has led to what some call "the x-ray cluster baryon 
crisis" [22] [23]. In fact, the gas in the central part (until a radius of about lh^ 1 
Mpc) of the cores of rich clusters is seen through the x-rays emitted by the gas 
in hydrostatic equilibrium and it has been recently seen that the fraction of the 
total mass in gas in those regions is large, f = {M gas /M to tai) — (0.05 — 0.10)/i~ 3 / 2 . 
Since the clusters in question are large clusters, one can think that this ratio is 
representative of the baryon fraction in the Universe Pb/pm = ^b/^dm- If so, 
Hdm = ^b/^ 1 and using the measured large value of f and the BBNS upper bound 
of 0.02/i" 2 on Q B one gets an upper limit of Q < (0.2 - OA)^ 1 / 2 [22]. If h > 0.16 
(as all measurements confirm, see section 3.4) and the bound on il B from BBNS is 
correct, this would mean that either we live in an open Universe (if A = 0) or in a 
Universe with A ^ (if we want it flat we need Qdm + = 1) or both. 



4. Thermal History Of The Universe 



As we have seen in sections 3.3 and 3.7, the present energy in radiation f2 7 ~ 1CT 5 , 
is much smaller than that of matter Q ~ ^dm — 0.1 — 1. This still holds true 
even when adding the contribution of relativistic neutrinos, with which Qr(T ) = 
4 x 1CT 5 /i~ 2 (g*/3.36), where g* are effective relativistic degrees of freedom and 
for photons and three relativistic neutrino species, g* = 3.36 (see section 5). Due 
to the different evolution with temperature of the density of matter and radiation, 
Pr ~ T 4 and pu ~ T 3 , the radiation was dominant in the past at T > T eq , 
where T eq is the temperature at matter-radiation equality, p R (T eq ) = pM(T eq ), 
T eq ~ 5.8eV Q H 2 (3.36/g*(T eq )) [5]. 

At even earlier times, particles of mass m could be produced and formed part of 
the "primordial soup" together with photons when T > m. When the temperature, 
that is the characteristic energy of the photons, was larger than the binding energy 
of a certain system, this system could not survive as such. This happened with 
atoms for T > 1 eV and with nuclei for T > 1 MeV. This mean that atoms became 
stable for the first time at the recombination epoch, T ~ 1 eV and nuclei were first 
formed at the BBNS epoch, T ~ 10 - 0.1 MeV. 

Going to still earlier times and larger T, QCD predicts that chiral symmetry 
breaking and the confinement of quarks within hadrons should happen at T ~ 100 
MeV. Even earlier, the electroweak symmetry should be restored at T > 100 GeV. 
Still earlier, we would encounter physics beyond the standard electroweak model, 
whatever it might be, until the Planck scale T ~ 10 19 GeV above which a quantum 
theory of gravity would be required. 

5. Relic Abundances of Particles 

Let us follow the histories of particles that were in equilibrium at a high tem- 
perature T in the primitive Universe. Stable particles X in equilibrium have the 
following number densities nx with respect to photons n 7 : 

—nx/n-y = /x(<7x/2), for relativistic particles, (mx «T), where fx — 1 for a 
boson and 3/4 for a fermion, and gx is the multiplicity of spin states (gx = 1 for a 
real spin boson, gx = 2 for the photon and a Weyl or Majorana fermion, gx = 3 
for a massive gauge boson, g x = 4 for a Dirac fermion); 

-n x /n y = (ir/8) 1/2 (gx/2((3))(m x /T) 3/2 exp(-m x /T), for non-relativistic parti- 
cles (m x » T). 

Particles go out of equilibrium when their rate of interaction V = nav becomes 
smaller than the rate at which T decreases, that is the rate of expansion of the 
Universe, the Hubble constant H, since —T/T = a/a = H (or, equivalently, when 
the mean interaction time T~ l becomes larger than the age of the Universe H^ 1 ). 
Particles for which the condition T > H was never fulfilled were never in equilibrium 
and their number has to be computed in other ways (this is the case of axions, for 



example). Different particles go out of equilibrium at different "freeze-out" or 
"decoupling" temperatures T/. D ., depending on their interactions. 

Relativistic particles with weak interactions go out of equilibrium in the early 
Universe at 7/. . ~ 1 MeV. The estimate of this temperature is very easy. For 
light left-handed neutrinos vl, for temperatures smaller than the mass of the vector 
bosons that mediate the interaction, the W and Z, (M v ~ 100 GeV) but larger than 
any other masses (of neutrinos, leptons, quarks ...) the cross section on dimensional 
grounds is a ~ G 2 F T 2 , with G F the Fermi constant, G F ~ 10~ 5 /GeV 2 ~ My 2 . The 
number density of light particles is n ~ T 3 , and their velocity v ~ c = 1. On the 
other hand, H ~ \[G T 2 , with G the gravitational constant, G ~ Mp 2 . Thus the 
interaction rate goes as T ~ nav ~ G F T 5 , it decreases faster with T than H and 
at T = T/.o. T becomes smaller than H. By equating T = H the mentioned T/. is 
obtained. 

The number n„ of light neutrinos (that are relativistic at T ~ lMeV) per co- 
moving volume is kept constant after their freeze-out, but the number of photons 
is increased by the annihilation of e + e~. When T drops below a mass threshold 
for the production of pairs of certain particles, e + e~ in this case, these particles 
cease to be produced and only annihilate "heating" the photons and all the other 
interacting particles but not the decoupled species. Consequently, the tempera- 
ture of the decoupled particles, T v in this case, becomes smaller than the photon 
temperature T. Actually, it is not that the temperature of photons increases, it 
just decreases at a smaller rate for a while. We can compute the ratio T u /T using 
entropy conservation. In fact, in thermal equilibrium the total entropy per comov- 
ing volume, S = sa 3 , is conserved. Here s is the entropy density s = {p + p)/T, 
thus s = (4/3)(p/T) for radiation (and relativistic particles). Let us compute the 
entropy before and after e + e~ annihilation and use entropy conservation, with the 
approximation that the scale factor a before and after is the same, namely as- 
suming there is not enough time for a to change much. The entropy before the 
annihilation has a contribution due to e + e~, s e = 4 x (7/8) x (27r 2 /45)T 3 and a 
contribution due to the photons s 7 = 2 x (27r 2 /45)T 3 . After the annihilation only 
the photons remain with T = T aft cr > T bcforc , since (2 + 7/2)Tg eiore = 2T 3 fter . The 
temperature of the neutrinos is not changed (because neutrinos are decoupled) and 
remains equal to T before , T v = T beforc = (4/11) x / 3 T. Thus, the number of neutri- 
nos plus anti- neutrinos, n^ u+p ^ = rii of every light species z/ e , z/ M , u T , that was before 
rii/n^ = (3/4:)(g u /2) becomes nj/n 7 = (3/ll)(^/2), that is, n, = (^/2)(115 /cm 3 ), 
because the number of neutrinos even after their decoupling is still proportional to 
T 3 (see e.g. [5]). 

Knowing T v we can compute the contribution of each z/-species that is still rel- 
ativistic to the present radiation energy density pn, that is usually parametrized 
as p R = (tt 2 /30) <?*(T) T 4 . Here g*(T) = Eb9b(T b /T)* + g F {T F /T) 4 is the 
effective number of relativistic degrees of freedom at the temperature T, the sums 
run over all relativistic bosons B and fermions F at T, gs and g F are the number 
of bosonic and fermionic degrees of freedom and Tg and T F are the temperatures 
of each species, which are equal to T only for species still in thermal equilibrium. 



For every v species g*(T ) = 2 x 7/8 x (T v jT f = 2 x 7/8(4/1 1) 4 / 3 = 0.454. 

Assuming all three neutrino species are still relativistic, photons and neutrinos 
account for g,(T ) = 2 + 3 x 0.454 = 3.36, and Q R h 2 = 4 x KT 5 (#*/3.36) (as 
already mentioned in section 4). However, it is possible that one or more neutrinos 
are non-relativistic at present (if m v > T Q = 2.3 x 10~ 4 eV), since the experimental 
upper bounds are m Ue < 5 eV, m Vji < 160 keV, m Vr < 24 MeV. In this case, the 
contribution of the non-relativistic light neutrinos to the density of the Universe 
now is p u = X^=i m ^ n i = ^uPc which means 

3 2 

£> Vi = 92eV n„h 2 (—) . (11) 

i=l 9v 

The sum runs over all non-relativistic neutrino species lighter than 1 MeV (because 
we use n v for neutrinos relativistic at freeze-out). 

From Eq. (11) an upper bound on neutrino masses follows from an upper bound 
on the present density of the Universe, f2„ < VL . The best bound on il h 2 comes 
from a lower bound on the lifetime of the Universe (see section 3.6). In a matter 
dominated Universe (with A = 0), t >, 1.3 x 10 10 yr requires Q Q h 2 < 0.4, but 
adding the reasonable bound h >^ 0.5 one obtains Q Q h 2 < 0.25. These two bounds 
imply (using Eq. (11)) m Vi ^ 37 eV and J2i < 23 eV, respectively [24]. With 
a non-zero cosmological constant, A > 0, these bounds are relaxed a bit, since Q h 2 
can be larger for the same t D . These bounds apply to the masses m Vi of neutrinos 
with two degrees of freedom, i.e. Weyl of Majorana, and full weak interactions, 
thus g v = 2. The right-handed neutrino components, even when they exist, are not 
counted, because either they are never in equilibrium or decoupled much earlier 
than the left handed components, and are therefore much less abundant. 

Heavy relics X of mass mx that become non-relativistic while still in equi- 
librium, have a different history. While annihilations are in equilibrium, i.e. 
^annih. > H, the number of non-relativistic particles decreases with the Boltzmann 
factor nx/T 3 ~ exp(—mx/T). When annihilations cease (because the rate T ann ih. 
becomes small with respect to H), nx/T 3 departs from its equilibrium value and 
soon it becomes constant. The larger the annihilation cross-section is, the longer 
annihilations remain in equilibrium, and the smaller is the relic abundance Q, thus 
[25] 

nh 2 * 1 x 10 ~ 3W . (12) 

<a a v> 

Here < a a v > is the thermal average of the non relativistic annihilation cross- 
section that is always proportional to i> _1 , and <> indicates the average over a 
thermal distribution of momenta (average over initial states and sum over all final 
states). 

Eq. (12) implies that particles of mass mx of order GeV- TeV with annihilation 
cross sections of "weak" order, i.e. < a anni hV >~ G 2 F m 2 x ~ 10~ 37 cm 2 (mx/GeV) 2 



for m x < M z and < cr annih v >~ a/m 2 x ~ 10~ 37 cm 2 (mx/TeV)~ 2 for m x » M z , 
have relic abundances Vt of order one, and may be, therefore, good DM candidates 
(a is the electromagnetic constant). 

Since elementary particle physicists are exploring extensions of the electroweak 
model, plenty of new particles with masses in the few GeV to a few TeV, are 
hypothesized. Many of them have interactions mediated by particles in the same 
mass range, what means cross sections of the weak order and, thus, they can account 
for n ~ 0.1 to 1 and are good DM candidates. This is the case of the Lightest 
Supersymmetric Particle, for example. 

In the above discussion no asymmetry between X-particle and antiparticle num- 
bers was assumed. A cosmic asymmetry would insure that annihilations stop when 
the minoritary population is depleted, leaving many more particles than expected 
without an asymmetry. 

6. Primordial Nucleosynthesis, Also in Crisis? 

When the temperature of the Universe became smaller than the binding energy 
of nuclei, at T ~ 0.1 MeV, nuclei first became stable. Because neutrons are heavier 
than protons n n /n p ~ exp[— (m n — m p )/T] < 1, there are less neutrons than pro- 
tons. Most of the neutrons end up into 4 He (approximately 25% of the mass of the 
Universe) together with an equal number of protons. The remaining protons stay 
as H (approximately 75% of the mass of the Universe), and trace amounts of D, 
3 He and 7 Li are produced (with n Dj 3 Hc /n H — O(10~ 4 ) and (n7 Li /n H ) — O(10~ 10 ). 
Primordial abundances of the light-elements are difficult to infer from observations 
and experts differ on the relative weight they give to different measurements and 
methods of inference. While already in 1964 it was known that stars can only 
produce less then 5% of the existing 4 He, only in the early 70's it was shown that 
D, being very weakly-bound, is easily destroyed in stars (mostly into 3 He) but it 
cannot be produced. 

The predictions of BBNS depend on the value of the baryon-to-photon ratio 
7] = nB/n 7 during the nucleosynthesis, and it is remarkable considering that the 
abundances are vastly different, that there is a range of rj ~ O(10~ 10 ) for which 
realistic abundances are obtained for all the light element. The prediction of the 4 He 
abundance also depends of the number of relativistic neutrino species in equilibrium 
during the BBNS, N v . The abundance of 4 He increases with N v . Thus present 
upper bounds on 4 He give an upper bound on N v of about 3. N v parametrizes 
also any non-standard contribution to the energy density of the Universe during 
the BBNS, so it would be easy to explain a BBNS prediction of N v larger than 3. 
However, a bound N v < 3 would require at least one of the three known neutrinos 
to contribute less than a relativistic neutrino during BBNS. For example, a heavy 
neutrino could decay before the BBNS, i.e. with a lifetime < 10 sec. Otherwise, one 
would suspect that there is something wrong in the arguments leading to the BBNS 
bound on N v . Precisely a bound N v < 2.6 (95% C.L.), with a value N u = 2.1± 



0.3 (1 a error), is the claim of a recent paper entitled "Big Bang Nucleosynthesis 
in Crisis" [26]. One of the main differences in this paper is a lower assumed range 
for D+ 3 He. Because the abundance of 4 He increases with 77 while that of D and 
3 He decrease with 77 ( 7 Li has a dip, but increases with 77 for 77 >^ 3 x 10 -10 ) a lower 
assumed range for D+ 3 He, pushes the acceptable range of 77 in this paper, to larger 
values. Fig. 3 (taken from [5]) shows the dependence of the different abundance on 
77 and of the 4 He abundance on N v . Because the abundance of 4 He also increases 
with N v , in order to prevent getting too much 4 He with the larger values of 77 
selected by the lower D abundance, N v must be smaller. Given the dispersion in 
different measurements of primordial abundances (26), to claim, a "crisis" seems 
premature. However one should keep in mind that the primordial abundances of 
4 He, D+ 3 He and 7 Li leave a narrow range of 77 where agreement with predictions 
for all of them is achieved and when observational uncertainties decrease further in 
the near future (mainly in the primordial abundance of D [27]) this narrow range 
will be further squeezed. 

Let us briefly see where the different dependences with 77 and N v of the various 
abundances shown in Fig. 3 come from (for more details see e.g. [5]). Let us start 
our description at t ~ 1CT 2 sec, when T ~ 10 MeV, 7, v, P, e + , e~, n and p 
are in equilibrium and n n /n p ~ 0.9 due to the difference in mass between protons 
and neutrons m n — m p = 1.293 MeV. At about t ~ 1 sec, T ~ 1 MeV, weak 
scatterings and annihilations go out of equilibrium, thus z/'s decouple and the ratio 
n n /n p freezes-out at n n /n p ~ 1/6 and decreases from then on only because neutrons 
decay (thus the BBNS predictions depend on the neutron lifetime). 

The exact temperature at which the weak decoupling occurs depends on the 
value of Hubble parameter H. The larger H, the sooner the decreasing rate of 



FIGURE 3. Predicted primordial abundances as function of rj and iV, 



weak scatterings and annihilations becomes equal to H, thus the weak decoupling 
(or freeze-out) occurs earlier, consequently the n n /n p ratio at freeze-out is larger 
and, because practically all neutrons available when nucleosynthesis finally happens 
end up in 4 He, the predicted abundance of 4 He is also larger. Since the density 
of the Universe p increases with the number of neutrino species N u , H ~ ^fp and 
consequently the predicted abundance of 4 He also increase with increasing A^. 

Even if the temperature at this point, T ~ 1 MeV, is lower than the binding 
energy of 4 He, the synthesis of this nucleus (and consequently, that of heavier 
nuclei) cannot happen because the deuterium D has a very small binding energy, 
i.e. it is constantly destroyed by photons with E > 2.2 MeV. Because there are so 
many more photons than baryons, we need to wait until the amount of photons in 
the tail of the energy spectrum with E > 2.2 MeV becomes smaller than the number 
of D nuclei, namely when n 7 exp(— 2.2 MeV/T) < When this happens, at T ~ 
2.2MeV/^n(77 -1 ) the so called "D bottleneck" finishes and 4 He is formed. Because 
until the end of the "D-bottleneck" period neutrons continue decaying, a large value 
of r\ leads to a shorter time for neutrons to decay what leads to larger amounts of 
4 He produced. In order to understand the duration of this period it is necessary to 
take into account the reheating of photons due to e + e" annihilation (see section 5) 
that happens precisely in this period. Due to this effect the temperature decreases 
less rapidly and the "D-bottleneck" finishes at t ~ 200 sec (without reheating it 
would finish at t ~ 7.35 sec instead). 

At this point practically all neutrons end up in 4 He thus n( 4 He) = n n /2. Using 
n n /n p ~ 1/7 we obtain with a back-of-an-envelope calculation the right amount of 
4 He, V( 4 He)=0.25 : 

4(n n /2) = 2(n n /n p ) 
[4n( 4 He) + n(H)] m p 4(n n /2) + n p - n n 1 + (n n /n p ) 

that is precisely 0.25. Here n( 4 He), n(H) are the number densities of 4 He and H 
nuclei and m p stands for the mass of both, protons and neutrons. 

Trace amounts of 7 Li are also produced, but not heavier elements. The elements 
with atomic weight 5 and 8 are not stable and the gap in mass is then too large for 
heavier elements to be formed. This gap is bridged in stars at higher T though three- 
body collisions. There are two different processes through which 7 Li is produced, 
that cross in a dip at rj ~ 3 x 10~ 10 . 

In addition, some amount of D and 3 He are left unburnt, as the reactions that 
burn them into 4 He freeze-out. The rates of these reactions are proportional to rj. 
So for larger values of rj the reactions that burn D and 3 He are longer in equilibrium 
and less D and 3 He remain. 

A range in 77 is translated into a range of ilsh 2 , since 77 = (%/n 7 ) = 
(n B p c /m p )/(2£(3)T 3 /7r 2 ). With T = 2.73°K one obtains r/ = 272.2 x 10" w Q B h 2 . 
Bounds on 77 have changed over time, mainly due to observational uncertainties 
and different treatments of the errors. The ranges of the "Chicago group" have 
changed from r] = (4 - 7)10" 10 , i.e. 0.015 < Q B h 2 < 0.026 (and N u < 3.4) in 
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1980's, to rj = (2.8 - 4)l(r 10 , i.e. 0.012 < VL B h 2 < 0.015 (and N v < 3.3) in 
1991, to r) = (2 - 6.5)10~ 10 , i.e. 0.007 < Q B h 2 < 0.024 (and N v < 3.9 at 95% 
C.L.) in 1995 [28]. The newer "Ohio group" prefers r] = (3.8 - 5.2) lO" 10 , i.e. 

0. 014 < Q B h 2 < 0.019 (and N v < 2.6 at 95% C.L. as mentioned above) [26]. 

As a final comment, let us mention a particularly intriguing explanation for the 
result N u < 3, if confirmed. The solution may be an unstable v T with mass in 
the MeV range decaying with a lifetime r < 10 sec into invisible particles [29], 

1. e. z/'s or i/(f> with <fi a Goldstone boson, a Majoron. In particular if the daughter 
particles include u e and m UT = 1 — 10 MeV, these extra v e would deplete n n due 
to the process u e n — > p e~, because the opposite reaction V e p — > n e + requires an 
energy of the daughter neutrinos E De = m UT /2 > m n — m p ~ MeV and is suppressed 
for m VT < 10 MeV. This effect can be played off against the increase in n n (and 
consequent increase in 4 He) due to extra contributions to the density p during the 
BBNS to allow even for N v = 16 ! 

7. Structure Formation 

The Universe looks lumpy at scales A ~ 100 Mpc, we see galaxies, clusters, super- 
clusters, voids, walls. But it was very smooth at the surface of last scattering of the 
CMBR (i.e. at electromagnetic decoupling, when ions and electrons first formed 
atoms) and later (see section 3). Inhomogeneities have been seen as anisotropies 
in the CMBR, so (Sp/p) cannot be much larger than ST /T < 10~ 5 . So inho- 
mogeneities in density start small and grow through the Jeans (or gravitational) 
instability: gravitation tends to further empty underdense regions and to further 
increase the density of overdense regions. One can follow analytically the evolution 
due to gravity of the density contrast 5p = (p(x) — p)j p (where p is the average 
density) in the linear regime where Sp/p < 1. In a static fluid the rate of growth 
of Sp/ p is exponential, but in the Universe (an expanding fluid) it slows down into 
either a power law, Sp/p ~ a, in a matter dominated Universe or it stops, Sp/p ~ 
constant, in a radiation or a curvature dominated Universe. We have not so far 
mentioned the possibility of a Universe where the curvature term dominates the 
r.h.s. of Eq. (6). Using this equation one can see that a matter dominated open 
Universe becomes curvature dominated for a > Q /(l — £l ) and no further growth 
of the contrast of perturbations in matter can occur. 

Perturbations have different physical linear dimensions A which increase with the 
Hubble expansion, A = a(t)X com . Here a(t) is the scale factor and A com is the linear 
dimension measured in comoving coordinates (those that expand with the Hubble 
flow). With the usual choice of a = 1 at present, A com is the present actual linear 
dimension. Since a ~ t a with a < 1 (see Eq. (8)) while the horizon ct grows linearly 
with t, the horizon increases with time even in comoving coordinates, encompassing 
more material as time goes. When A = ct we say the perturbation of size A "enters" 
into the horizon, we could better say that the perturbation is first encompassed by 
the horizon. This moment is called "horizon-crossing" and it happens at different 



times for different linear scales A, larger scales cross later. 

What is the origin of these primordial fluctuations? In the standard BB model 
perturbations can at most cross the horizon only once, coming from outside the hori- 
zon, where they cannot be generated because there can be no causal interactions. In 
inflationary models each perturbation crosses the horizon twice. Inhomogeneities in 
density are generated as quantum fluctuations, expand exponentially going outside 
the horizon during inflation and re-enter later, after inflation is over. Inhomo- 
geneities could also be generated within the horizon by defects, such as strings and 
textures, that are remnants from phase transitions in the early Universe. 

It is convenient to specify the spectrum of fluctuations at horizon-crossing, 
(5p/p)hor- A Harrison-Zel'dovich spectrum is scale invariant at horizon-crossing, 
namely (dp/ 'p)hor = constant. COBE observations have shown this spectrum is in 
fact scale invariant or very close to it. 

After horizon-crossing, physical interactions act upon the inhomogeneities and 
generate a "processed" spectrum. Three main processes compete. Gravitational 
interactions tend to produce a gravitational collapse, pressure changes tend to 
produce pressure waves, and there may be damping, caused by the free streaming 
of non-interacting particles (such as relativistic neutrinos) or due to collisions (this 
is the case of the Silk damping for baryonic matter, due to photon viscosity). A 
gravitational collapse only occurs if the distance a pressure wave travels before the 
collapse, f sound ^free-fall = Aj eans , is smaller than the size A of the inhomogeneity, i.e. 
A > Ajeans- If so, roughly speaking, pressure waves cannot scape from the collapsing 
inhomogeneity. This is Jeans criterium for a collapse. 

For baryons before recombination i> SOU nd — c/v^ (since v^ ound = dp/ dp and the 
pressure is due to photons, thus p = p/3), thus the Jeans length is larger than the 
horizon and the density contrast (Sp/p)B of perturbations in baryons within the 
horizon cannot grow. After recombination, matter (atoms) decouple from radia- 
tion (since atoms are neutral), the pressure is due to non-relativistic H atoms thus 
the sound speed and consequently the Jeans length decrease steeply and pertur- 
bations in baryons can star growing. Moreover, Silk damping at decoupling erases 
inhomogeneities in baryons smaller than those corresponding to clusters of galaxies. 

We can now see that structure could not have been formed with baryons alone. 
Because the Universe has become matter dominated at t eq before recombination, 
t eq < t rec , perturbations in baryons can grow as (Sp/p)s ~ o, from recombination 
onwards. However, the maximum growth factor since then is a /a rcc pa 10 3 , and 
with (Sp/p) rcc < 1CT 4 , as shown by COBE and other CMBR anisotropy experi- 
ments, (5p/p)b at present could not have reached 1. However, collapsed objects 
(galaxies, clusters etc.) have Sp/p > 1. Therefore, we conclude that there is not 
enough time from recombination to the present for inhomogeneities in baryons to 
grow to (Sp/ ' p)b — 1 (after that, in the "non-linear" regime, Sp/p grows very fast) 
unless baryons after recombination (i.e. when finally they can collapse into struc- 
tures) fall into already existing potential wells. These potential wells can only be 
formed by matter not coupled to photons, namely dark matter (DM) . Perturbations 
in the DM can start growing before recombination, as soon as the Universe becomes 



matter dominated (by the DM itself), so they have enough time to become non 
linear by now. By recombination the DM has formed the potential wells into which 
baryons fall. The necessity of a matter component in which density perturbations 
at the moment of recombination may be much larger than the perturbations in 
baryons, is a strong motivation for the existence of DM. 

The question of which structures are formed first leads to the distinction of three 
types of DM: hot, warm and cold. The following naive picture of how fluctuations 
evolve helps understanding the differences of these three types. 

Let us define as t ga i, the moment at which the growing horizon encompasses for 
the first time a perturbation of the size of a typical galaxy that contains, 10 11 M Q 
(including the dark halo), where M is the mass of the Sun. The temperature 
of the Universe at that moment turns out to be T ~ 1 keV. Hot DM particles 
are still relativistic at t ga x, i.e. their mass is m < 1 keV . In this case, at t ga i 
each of the particles in a fluctuation of the size of the galaxy has moved from its 
original position a distance ct equal to the size of the fluctuation. The volume 
occupied by the particles will expand with the horizon. Thus, a moment later 
the fluctuation is erased. This is damping by free streaming. At the moment 
when the DM particles become non-relativistic, their motion becomes negligible. 
Thus the smallest structure that can survive the "free streaming" of relativistic 
particles and grow, is that one encompassed by the horizon when the particles are 
becoming non-relativistic. Light neutrinos are hot DM and the smallest structure 
that can survive contains the Jeans mass, Mj eans = 3 x 1O 15 M / \m u /30eV) 2 that 
corresponds to a supercluster. Therefore, with hot DM (HDM) (such as light 
neutrinos of m < 30eU) superclusters of galaxies form first and galaxies must form 
later, through the fragmentation of the larger structures. Simulations have shown, 
however, that there is not enough time to form galaxies. Thus HDM is rejected as 
the dominant DM component. Particles of mass m ~ keV, that are just becoming 
non-relativistic at t ga \, are warm DM. With warm DM galaxies form first, but 
barely. Cold DM (CDM) is such that perturbations smaller than a galaxy, even of 
the size of the smallest dwarf galaxy (1O 6 M ), are not erased and can grow. 

Simulations have shown that CDM must be the most abundant form of matter, 
because the "processed" spectrum of perturbations (i.e. the present spectrum once 
the effect of physics acting within the horizon has been considered) generated in 
standard CDM models reproduces the observations within 10%. Standard CDM 
models, make the simplest assumptions, namely I^cdm + — f2 = 1, A = 0, scale 
invariant perturbations at horizon crossing, and a scale independent "biasing" by 
which only the highest peaks in the CDM density distribution end up forming 
galaxies. There is only one feature in the processed spectrum of CDM perturba- 
tions, a change of slope at the present scale that corresponds to the horizon at the 
moment of matter-radiation equality, A eq . Perturbations with A < A eq entered into 
the horizon at t < t cq , when the Universe is radiation dominated. They cannot 
grow while the Universe is radiation dominated, so they all start growing together 
at t = t cq and they roughly have the same amplitude today if they all started with 
the same amplitude at horizon crossing. Perturbations with A > A cq , instead, en- 



tered into the horizon at t > t cq , when the Universe is matter dominated, so they 
started growing immediately. Consequently, perturbations at larger scales entered 
later, had less time to grow and their amplitude is smaller at present. Once A eq (the 
location of the change slope) is fixed, the only remaining free parameter is an over- 
all normalization of the CDM spectrum, now provided by the CMBR anisotropy 
measured by COBE at large scales, 9 > 20°. Density perturbations at these large 
scales entered into the horizon very recently (so they did not grow much), thus 
providing a measurement of {Sp/p) at horizon crossing, (Sp/p)h OT . For more details 
on this section so far see e.g. [5]. 

While both the shape and normalization so obtained are almost right, they do not 
fit the observations [30]. In Fig. 4 [31] the power spectrum of density perturbations 
predicted by standard CDM models (solid line) normalized by COBE data (the box 
on the l.h.s.) and other data (the points with vertical error bars) are shown. As 
can be seen in the Fig. 4 (solid line) the spectrum of standard CDM models has 
too much power on small scales (large k~ A^ 1 ), the scales of galaxy clusters and 
smaller. Also in Fig. 4 one can see the failure of HDM to account for the data (see 
the short-dashed curve) mentioned above. 

Once the normalization given by COBE is fixed, there are several possibilities 
to change the spectrum to agree with observations. Because HDM tends to erase 
structure at small scales, one of the solutions consists in adding to the CDM a bit 
of HDM, namely neutrinos, in what are called mixed DM (MDM) or hot-cold DM 
(HCDM) models [32]. See the dotted line in Fig. 4. In particular, models with 
Q u = 0.2, what corresponds to J2i m ui = 5 eV, and the rest of fl Q in CDM plus 
some baryons, with £1 = 1 and A = works well. Another possibility is that of a 
"tilted" primordial spectrum of fluctuations at horizon crossing, one that slightly 
favors larger scales over smaller scales (instead of the flat, scale invariant, Harrison- 
Zel'dovich spectrum) within the COBE observational limits. This is called "tilted" 



FIGURE 4. Comparison of the measured power spectrum of density perturbations and the 
predictions of several DM models (notice that A -1 ). 



CDM (TCDM) [33]. See the long-dashed line (labelled with n=0.8) in Fig. 4. 

Another family of solutions is obtained by realizing that a shift towards larger 
scales of the only feature in the CDM spectrum, i.e. A eq , the scale where the slope 
changes, is enough to provide good agreement with observations, since it effectively 
amounts to increasing the power of the spectrum at scales larger than the break 
point (A > A cq ) with respect to those smaller than it (A < A eq ). The break in 
the slope occurs at the moment of radiation- matter equality, p ra d = Pmattcr, when 
the Universe becomes matter dominated (by the DM). The scale of the break, A cq , 
corresponds to the present size of the horizon at matter radiation equality, ct eq . It 
is A cq ~ lOMpc (g„,/3.36) 1//2 {Vt h 2 )~ l , where Vt Q is essentially the present matter 
density, since Qr << Q q (recall that the total density is Q = Q Q + fl vac , if the 
cosmological constant A ^ 0) and is the effective number of relativistic degrees 
of freedom at t cq , g* = 3.36 for photons plus three relativistic neutrino species (see 
section 5). 

The equation A cq = (lOh^ 1 MpcjT^ 1 defines the "shape parameter" [34] T = 
^o^(5 , */3.36) -1 / 2 (notice that distances determined through red-shifts are given in 
h^Mpc units (see Eq. (5)). The data require T ~ 0.25±0.05, while standard CDM 
models (with the standard choices of h = 0.5, Q Q — 1, g* — 3.36), has T = 0.5. In 
fact, as we have explained a larger A eq , thus a smaller T, would provide agreement 
with data. In order to lower the value of T with respect to that of the standard 
CDM model one needs to either: 1- lower h [35] (what implies an older Universe), 
or 2- increase g* (namely increase the radiation content of the Universe at t eq ), or 
3- lower Q (i.e. take Q Q < 1), so that we either live in an open Universe (open 
CDM models, OCDM) if A = 0, or in a Universe with a cosmological constant that 
provides f2 vac = 1 — Q Q (ACDM models [36]), for example with h—1, ^cdm = 0.18, 
il B = 0.02, f2 vac = 0.80, or 4- a combination of all three above. We have already 
examine the consequences of the 3rd. solution in relation with the possible "age 
crisis", in section 3.6. Let us examine the consequences of the other two. 

If we want to keep Q Q = 1, with the standard value of g* = 3.36 one would 
need h = 0.3 to get T = 0.3, but this would lead to a very old Universe with 
t = 2.2 x 10 10 yr. If instead we keep the standard value of h i.e. h = 0.5 (this 
low-value of h is necessary with Q Q = 1 to account for t = 1.3 x 10 10 yr), the 
value g* = 9.33 is necessary to get T = 0.3, what amounts to the equivalent of 16 
relativistic neutrino species. However, primordial nucleosynthesis does not allow 
for much more than 3 neutrino species (except with a heavy unstable z/ T , see section 
6 and below). 

A way to obtain the large amount of radiation needed is through a heavy neutrino 
decaying into radiation, with the right combination of mass and lifetime, in so-called 
rCDM models [37] [38]. A massive neutrino matter dominates the energy density 
of the Universe as soon as it becomes non-relativistic m v >^ T (since n u ~ n 7 and 
p u = n u m u , p rad ~ n 7 T), thus their decay products radiation-dominate the Universe 
at decay. For m v < lMeV the right mass-lifetime combination lie on a narrow strip 
around the (usually called) "galaxy formation" bound. This bound [39] is obtained 



by requiring that the decay products of a massive neutrino cease to dominate the 
energy density of the Universe before baryons should start falling into gravitational 
wells, at recombination t rec ~ 10~ 5 t o , what yields r < (92eV/m^) 2 (r2 o /i 2 ) 2 10~ 5 t o . 
Near this bound, at the boundary between being irrelevant and harmful, unsta- 
ble neutrinos could help in the formation of structure in the Universe. A heavier 
neutrino, of m v ~ 1 — 10 MeV, necessarily v T) decaying at or just before nucleosyn- 
thesis, r = 0.1 — 100 sec, would also provide a solution [38]. If r— 0.1-few sec and 
v e are among the decay products, neutrons would be depleted, leaving room for the 
presence of up to 16 equivalent relativistic neutrino species during nucleosynthesis 
[29], as explained at the end of section 6. This v T decays too fast to get to matter 
dominate so the decay products never dominate the energy density of the Universe 
and additional unknown radiation has to be added ad-hoc to increase g* as needed. 
If instead t v ~ 10 — 100 sec, the v T can dominate the energy density of the Universe 
before decaying and its decay products provide the necessary extra radiation. The 
v T decay modes involved here should all be into neutral particles, v T — > 3 v's or 
v T — > u(j), with a Majoron (a zero mass Goldstone boson) for example. All visible 
modes, i.e. producing electrons or photons, are forbidden in the necessary range. 

All these modified CDM models seem to be able to fit present data, however they 
predict very different patterns of acoustic peaks in the CMBR anisotropy power 
spectrum £(£+l)Ce (see section 3.3 and Eq. (10)). The location and relative height 
of the peaks between £ ~ 200 and several 1000 characterize different models. The 
detection of the CMBR anisotropies so far obtained with the COBE satellite and 
balloon experiments, does not discriminate between models, but the next generation 
of observations will do it. There are two approved new satellite experiment with 
sub-degree resolution, MAP in the U.S. and COBRAS /SAMBA in Europe. They 
will give results as shown in the Fig. 2 [13] together with the predictions of four 
different models, three of them mentioned above SCDM (standard CDM), ACDM 
(CDM with A > 0) and OCDM (open CDM), and one labelled with "strings", 
that corresponds to structure formation due to inhomogeneities caused in a phase 
transition that produced "cosmic strings" (vortex- like solitonic remnants). The 
dark band around the SCDM line shows the expected experimental error in the 
C/s for an experiment like COBRAS/SAMBA, with 10' beam-size and 10 fiK 
noise (that would allow to resolve anisotropies up to ~ 0.5° for fi = 1) [13]. It 
is clear from Fig. 2 that the different theoretical models shown could be easily 
distinguished. 

8. Detection of the Dark Matter in the Halo of Our Galaxy 

Our own galaxy is spiral. The Sun is part of a thin disk of about 12 kpc of 
radius consisting mainly of stars and gas clouds orbiting the galactic center. The 
Sun is at r Q ~ 8.5 kpc from the center and moving at about v c ~ 220 km/sec. 
The rotation curve of the Milky Way is flat at about v c out to as far as it is mea- 
sured, showing that our galaxy has a spherical or ellipsoidal dark halo. Assuming a 



spherical and isotropic velocity distribution, the density profile of this halo is usu- 
ally parametrized as p(r) = p D [(a 2 + r 2 )/(a 2 + r 2 )} where p Q ~ 0.01M o /pc 3 ~ 0.4 
GeV/cm 3 is the halo density in the solar neighbourhood and a is the core radius of 
the halo (p ~ const at r < a), that has been estimated at a ~ 2 — 10 kpc. There is 
also a central luminous component of the galaxy, the spheroid or bulge. The distri- 
bution and motion of the luminous matter trace the combined gravitational effects 
of the three components, the disk, the bulge and the halo. The halo dominates at 
large radius. 

In the last ten years, attempts have been made to detect the galactic DM in exper- 
iments that could reveal its nature. Because the DM in galactic haloes amount to 
^haloes — 0.02/i -1 (see section 3.7) nucleosynthesis bounds on the baryonic density 
Qb (see section 6) do not preclude the possibility that "dark" baryons constitute 
some or all the DM in galaxy haloes like ours. Candidates for dark baryons could 
be macroscopic objects up to the size of "brown dwarfs". These are stars too small, 
m < 0.08M Q , for their centers to reach a sufficiently high temperature to burn hy- 
drogen (as in a normal star). These candidates are called MACHOs (for massive 
compact halo objects) as opposed to WIMPs (for weakly interacting massive par- 
ticles), that are elementary particle candidates for the DM. Elementary particle 
candidates include also massive neutrinos (actually the only DM particle candi- 
dates known to exist, but not known to be massive as needed). The three known 
neutrinos, if stable, are hot DM candidates, and the energy they would carry is too 
small to allow their detection in the dark halo of our galaxy (but the measurement 
of their masses in the range 1-few 10 eV could tell). The theoretically preferred 
cold DM candidates are WIMPs, with masses in the range 1 GeV- 1 TeV, and the 
axion, a very light boson that appears in the Peccei-Quinn solution to the strong 
CP problem of elementary particles, with a mass m ~ 10 _5 eV. 

WIMPs appear in many extensions of the standard model of elementary particles, 
such as the minimal supersymmetric extension of the standard model, where the 
usual WIMP candidate is the lightest "neutralino" (the lightest fermionic partner 
of the neutral gauge and Higgs bosons). As we have seen (see Eq. (12)), it is 
enough for a particle in the GeV- TeV mass range to have annihilation cross sections 
of weak order in order to have a relic abundance f2 ~ 0(1). WIMPs are being 
search for in direct and indirect DM searches and accelerator experiments. If these 
particles compose the halo of our galaxy they may be detected directly, through 
the energy they may deposit in collisions with nuclei within detectors, or indirectly 
through their annihilation products, coming from annihilations in the Sun or Earth, 
where they may accumulate, or in the halo itself. Direct WIMP searches look for 
energies of order mv 2 ~ 0(10 keV) deposited in ionization and phonons, in kg- 
size cryogenic crystals, or producing scintillation or excitations of superfluids and 
superconductors. Indirect WIMP detectors search mainly for high energy neutrinos 
produced in WIMP annihilations in the Sun and the Earth, in large water or ice 
Cerenkov light detectors (such as Superkamiokande, AMANDA, DUMAND and 
NESTOR) or underground detectors (such as MACRO). For a review of particle 
DM searches see e.g. [40]. 



Axions could be detected through their coupling to two photons, looking for 
the resonant conversion of halo axions to microwave photons in a cavity with a 
magnetic field. The first detector of this type capable to actually test halo axions 
is based at Livermore [41] and has started taking data recently. 

The most fruitful DM searches so far, have been those looking for MACHOs 
though the gravitational microlensing [42] of background stars in the Large Mag- 
ellanic Cloud (LMC), an irregular galaxy satellite of the Milky Way. In a gravita- 
tional lensing event an intervening object located in the line of sight of a luminous 
object deflects the light coming from it, producing multiple images. In a microlens- 
ing event, these multiple images are superposed, what is seen as an increase in the 
brightness of the lensed object. Three independent groups the American-Australian 
MACHO, the French EROS and the Polish- American OGLE collaborations (that 
seem to be competing for the most politically incorrect name) have seen microlens- 
ing events. They were looking for (and found in 1993 [43]) events in which a star 
from the LMC or the bulge of our galaxy symmetrically brightens and fades, as a 
MACHO (or a not seen faint star) passes near its line of sight (one could describe 
the event as a symmetric "anti-eclipse") during a time At that depends on the 
lensing mass M, At = 14Odays(M/M ) 1/2 . So far about 90 events have b een seen 
in the direction of the central bulge of our galaxy (where most lenses are expected 
to be faint stars in the bulge itself) and about 10 towards the LMC. These numbers 
are about a factor of 5 too small from what was expected from the LMC if the dark 
halo of the Milky Way consisted of MACHOs and about a factor of 3 too large from 
what was expected from the bulge (see e.g. [44]). 

The number of events towards the LMC is too small to be compatible with a dark 
halo dominated by MACHOs (within the mass range tested M ~ 10~ 7 — 1M ). 
Machos could account for up to ~ 30% of the expected halo mass, unless the halo is 
much lighter than presently believed [45]. The large number of microlensing events 
towards the bulge, where the lensing objects are expected to be mostly faint stars, 
is much larger than expected on the basis of assuming a spherical bulge in the center 
of our galaxy, as seen in regular disk galaxies. The excess indicates that the bulge 
of our galaxy is a bar with its axis oriented towards the Sun [46], in agreement with 
other observations, what would yield a larger star density along our line of sight 
towards the bar. An alternative explanation could be that the mass of the disk is 
much larger than in standard galactic models. More data will allow to distinguish 
these possibilities. 
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